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I. INTRODUCTION 



A time-independent classical oscillator is exactly solved when the evolution of the po- 
sition and momentum in the phase space is completely known in terms of the initial data 
or constants of motion. A constant of motion is a function of the initial data in phase 
space, and the evolution of an oscillator is completely determined by any two independent 
constants. For any linear second-order system, such as a time-dependent classical oscillator, 
two invariants (constants of motion) are readily constructed from two independent classical 
solutions. The time-dependent quantum oscillator can be exactly solved in terms of the clas- 
sical solutions, as the invariant operator, which was first introduced by Lewis and Riesenfeld 
and is quadratic in position and momentum, is known explicitly in terms of the solutions of 
the corresponding classical oscillator [|TJ . The invariant operator generates the Fock space of 
exact quantum states of the Schrodinger equation. Recently the time-dependent quantum 
oscillator has attracted much attention with the discovery of the geometric phase by Berry 
to study the nonadiabatic geometric phase for various quantum states, such as Gaussian, 
number, squeezed, or coherent states, which can be found exactly 

In spite of intensive study of time- dependent oscillators, the issue of action-phase (angle) 
variables remains not fully exploited classically and quantum mechanically. In particular, 
the quantum action-phase operators are still elusive and not completely resolved since it 
was raised by Dirac ||]. For a time-independent classical oscillator the Hamiltonian itself, 
being a constant of motion, turns out to be an action variable, and its conjugate is the phase 
(angle) variable. The transformation from the phase space to the action-phase space is a 
canonical transformation. A generating function that yields the action-phase variables can 
be found even for the time-dependent classical oscillator || . On the other hand, for a time- 
independent quantum oscillator many different methods and schemes have been introduced 
to define the action-phase operators ||. Recently a candidate for a phase operator conjugate 
to the Hamiltonian (action) operator has been proposed for the time-independent quantum 
oscillator |7|]. 

The aim of this paper is to find the canonical transformations from the phase space to 
the invariant spaces and to define the action-phase variables for a time- dependent quadratic 
oscillator. For that purpose we introduce pairs of classical invariants that are linear in po- 
sition and momentum and are expressed in terms of the classical solutions, which realize 
explicitly the canonical transformation from the phase space to the invariant spaces. We find 
the most general three-parameter quadratic invariants formed from the linear ones, which 
exhaust all quadratic invariants and include the one by Lewis and Riesenfeld. In each invari- 
ant space we define the action variable, quadratic in position and momentum, and find its 
conjugate phase variable. Analogously by using pairs of linear invariant operators that con- 
stitute the annihilation and creation operators for the Fock space of exact quantum states 
HH, we find the action operator for the time-dependent quantum oscillator. The three- 
parameter invariant operator is transformed into the canonical form of a number operator 
through a unitary (Bogoliubov) transformation. Similarly the invariant operator obtained 
from one classical solution is transformed into the one obtained from another classical solu- 
tion through a unitary transformation. A number eigenstate of one pair is a squeezed state 
of the corresponding one of another pair. We extend the current scheme for defining the 
phase operator to the time-dependent quantum oscillator. 
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The organization of the paper is as follows. In Sec. II we find pairs of linear invariants for 
the time-dependent classical oscillator and show the canonical transformation from the phase 
space to the invariant space. In each invariant space we define the action-phase variables. 
In Sec. Ill we find the action operator and define the phase operator conjugate to the action 
operator. 



II. CLASSICAL ACTION-PHASE VARIABLES 



We consider a time-dependent classical oscillator which is described by the Hamiltonian 

(1) 



H{t) = —^p 2 + -^(pq + qp) + -4V, 



where X(t), Y(t) and Z(t) depend explicitly on time. An invariant for the Hamiltonian ([!]) 
is a constant of motion and obeys the equation 



j t I(t) = p(t) + {I(t),H(t)} PB = 0, 



(2) 



where { , }pb denotes a Poisson bracket. To find any invariant it is necessary to find a pair of 
invariants, say, I\ and I 2 , which are linear in position and momentum and are independent 
each other, because the invariants satisfying Eq. (0) satisfy the multiplication property 



|-(/lJ 2 )+{/l/ 2 ,#}pB 

at 



h + h 



d 

—I 2 + {I 2 , H} PB 
at 



0. 



(3) 



Hence any analytical function of the invariants I\ and I 2 is also an invariant. For that 
purpose, we introduce two invariants (see Appendix A) 



a x {t) = Ui(t)p - —[ui(t) - Yu x {t)]q, 
a 2 (t) = u 2 (t)p - -y[u2(t) - Yu 2 {t)]q, 
where u\ and u 2 are two independent real solutions to the classical equation of motion 



d ( u 
di\X / 



+ 



XZ -Y 2 + 



XY -XY 



X 



u 



0. 



(4) 



(5) 



where overdots denote d/dt. We require that these solutions are normalized to satisfy the 
Wronskian condition 



— (UiU 2 - UlU 2 ) = I. 

Then a\ and a 2 are canonical conjugates of each other, satisfying 

{a 1 (t),a 2 {t)}pB = 1. 



(6) 



(7) 



It should be remarked that there are as many pairs of independent invariants as pairs 
of independent solutions. These pairs of invariants constitute canonical conjugate pairs, 
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satisfying Eq. (|7|), if Eq. (|6]) is satisfied. As any two independent solutions to Eq. (|5|), say v\ 
and v 2, can be expressed in terms of the original pair, u\ and u 2 , these pairs are related to each 
other through a canonical transformation. In quantum theory the canonical transformation 
corresponds to a Bogoliubov transformation and has the meaning of squeezing the quantum 
states, which will be discussed Sec. III. 

For a time-independent oscillator with constant X, Y and Z, the Hamiltonian (|l|) is an 
action variable. For the time-dependent oscillator of this paper, the Hamiltonian cannot 
become an action variable, since dH(t)/dt ^ 0. Instead we look for the most general 
quadratic invariant of the form 

X = l -(B + A* + A)a\ + ^(A* - A)2 ai a 2 + l -(B -A*- A)a 2 2 , (8) 

where A is a complex constant and B is a real constant. In fact, Eq. (g) describes real 
three-parameter invariants that exhaust all real invariants that are quadratic in position 
and momentum. Through a canonical transformation, which is a rotation plus a scaling in 
the (01,02) space, Eq. ([Bp can be transformed into the canonical form 

T{t) = \[a\{t)+al{t)\. (9) 

The canonical transformation has the physical meaning of squeezing in quantum theory, 
which will be discussed in Sec. III. Constant values of the action (Q) depict circles in the 
(01,02) plane with polar coordinates defined through 



a x {t) = yj2X(t) cos 6{t), a 2 = yj2I(t)sm0(t). (10) 
The phase variable, simply given by 

a 2 (t) 



6(t) = tan" 1 



(11) 



Oi(t)_ 

up to arbitrary additive constant, is an invariant, and is the conjugate of the action (§): 

{l(t),6(t)} PB = l. (12) 

It is then trivial to show that 

{X(t),tan^(t)} PB = sec 2 ^(t). (13) 

The transformations from (q,p) to (oi,a 2 ) and from (01,02) to (X,9(t)) are canonical, so 
their volume elements are preserved: 

dq A dp = da\ A da 2 = dX A d9. (14) 

The volume of the phase space with X < X is then 

/ dqAdp=[ ddi A do 2 = / dX A d6 = 2ttX . (15) 

Jl<l Jl<l Jl<l 
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To make manifest the correspondence between the classical and quantum theory in Sec. 
Ill, we introduce another pair of invariants 



a{t) = {i)^u*(t)p- 
a\t) = {-i)iu{t)p- 



1 



X(t) 
1 

W) 



u*(t)-Y(t)u*(t)}q^. 
u(t)-Y(t)u(t)]q\, 



(16) 



where u is a complex solution to Eq. (||), normalized to satisfy the Wronskian condition 

1 , . 



— \uu — u u) 



(17) 



a and a* are the complex conjugates of each other and satisfy the Poisson bracket 

{a(t),a*{t)} PB = -i. (18) 

The factor — i is a consequence of complex nature of a and a*. If u{t) = [iu\(t) + u 2 (t)\/y2, 
the complex invariants (|I~6| ) are given by 



a 1 {t)+ia 2 (t) a 1 (t)-ia 2 (t) 
a{t) = -j= , a (t) = - 



V2 



V2 



(19) 



Hence there is a canonical transformation from the phase space (q, p) to the invariant space 
(a, a*). The action variable (0) is now given by 



T(t) = a*(t)a(t) 



u u\p + 



F 1 d 
X ~ Xdt 



ln(w*n 



,1/2 



~» 2 2 

9? + 



4u*u 



(20) 



A general quadratic invariant of a and a* can be transformed into the form ( p0[ ) through 
a canonical transformation, which has the physical interpretation of squeezing in quantum 
theory, as will discussed in Sec. III. 

As a* is the complex conjugate of a, the phase of a can be defined by 



a{t) = y/ltye"*®, a*{t) = ^l{t)e- i0a{t) . 
The polar form in Eq. ( pip leads to the phase variable in the form 

a {t)= l -[lna*(t)-\na(t)]. 

The equivalence between the phases (|TT]) and (|22|) can be shown using the identity 

~i + (a 2 (t)/ ai (t)y 



6(t) = tan" 1 



a 2 (t) 




ai(t). 
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i - (a 2 (t)/ai(t)) 



a (t). 



It follows that 9 a is the conjugate variable of I: 

{J(t) A(t)} PB = l. 



(21) 



(22) 



(23) 



(24) 
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After expressing q and p in ([Ll]) in terms of a and a* by using fll7|) , and writing the complex 
solution in a polar form 



u(t) =p(t)e* > »®, 



(25) 



we find 

q = ua + u*a* = 2\fTpcos(6 a — 6 U ) 



1 Y 2\fl 

P = -r7i.ua + u*a*) — — (ua + u*a*) = (p — Y p) cos(6> a — 6 U ) H sin(# a — 6 U ). (26) 

P 



X 



X 



X 



There is also another phase variable from the action (p0|). As Eq. (]20f) implies that 
constant X describes an ellipse in the (q,p) phase space, the phase variable can be defined 
by i 



cost? 



sin?? 



1 



Au*ul 



q, 



1 



p 



H_lAi n ( M * M )i/2 
X Xdt V ; 



Q ■ 



(27) 



The phase (|27|) is the conjugate variable of the invariant through the generating function 

" - 1- < 28 > 

where 

(29) 



S(T,q\t) = / p(I,q;t)dq. 

J qo ;X=fixed 



Hence the two phase variables (P7| ) and (|2l|) are related by 



(30) 



The canonical transformations from (q,p) to the invariant space of action-phase variables 
(X, 9) or (X, $) preserve the phase volume 



dq A dp = dl Ad6 = dl A d&. 



The area for X n < X is 



dqAdp= dl Add = dl Add = 2nl . 

I<I Jl<l Jl<l a 



(31) 



(32) 
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III. QUANTUM ACTION-PHASE OPERATORS 



In quantum theory invariant operators obey the quantum Liouville-von Neumann equa- 



tion 



ihp(t) + [i(t),H(t)] = o. 



(33) 



Corresponding to the pair of classical invariants (|T6|), there are a pair of invariant operators 



a+(t) = 



u*{t)p- 
u{t)p 



1 



X(t) 
1 



X(t) 



u*(t)-Y(t)u*(t)]q 
u(t) - Y(t)u(t)]q 



(34) 



where u is a complex solution satisfying (^]). Moreover, the invariant operators satisfy the 
standard commutation relation 



[a(t),a ] (t)] = 1. 



(35) 



There is the correspondence between the classical and quantum theory via a <-> a and 
a* <->• cr, except for some operator ordering problems. Just as the pair of classical invariants 
(p"6|) are fundamental in the sense that any analytic function F(a, a*) is an invariant, the 
pair of invariant operators (|34T) are fundamental, since any analytic function F(a, a*) is also 
an invariant operator. This fact is guaranteed by the multiplication property of solutions of 
the Liouville-von Neumann equation, 



d 











ih-(hh) + [hi 2 , H] = ih-hit) + [h, H]\i 2 + h[ ih—i 2 + [h, H 



dt 



dt 



dt 



0. 



(36) 



Hence the most general, quadratic, Hermitian, invariant operator is spanned by a) 2 {t), 
a*(t)a(t) +a(t)cv(t), and a 2 (t) with one complex and one real parameter, which corresponds 
to the classical one (|^). 

The most general, quadratic, Hermitian invariant that can be constructed from the pair 
(|34h has the form 



±^{t) + |[at(t)a(f) + a{t)a\t)) + ^-a 2 {t), 



(37) 



where A is a complex constant and B is a real constant. Another pair of invariant operators, 
say b(t) and b^(t), obtained by replacing u(t) by an independent complex solution to Eq. 
(|), say v(t), do not change the general form (p7|), since these are related to the original pair 
through the Bogoliuobov transformation 



b(t) = aa{t) + 0{t), 
&t(t) = a*a){t) +/3*a(t), 



(38) 



where a and (3 are determined by 
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v(t) = a*u{t) - 0*u*(t), 
v*(t) = au*(t) -f3u{t). 

A unitary (Bogoliubov) transformation can be found, 



(39) 



a(t) = cosh ra(t) + sinh re lS a(t) 

= S(r,8)ti(t)S\r,5), 
a^(t) = cosh ra (t) + sinh re~ lS a(t) 
= S(r,8)a(t)S\r,8), 

where the unitary operator is a squeezing operator 



l -{ z *h\t)-zh ] \t)} 



S{z) = exp 
which transforms (137T) to a canonical form 



z = re 



(5+71 



(40) 



(41) 



B t t 
At) = g- [a (t)a(t) +a(t)a\t)]. 

Here the squeezing parameter r, phase 5, and the coefficient B are determined by 

1 



(42) 



e i5 tanhr = — [-B ± V-B 2 - A* A], 

B = ^[A sinh 2re~ iS + B cosh 2r + A* sinh 2re i(5 ] . 



(43) 



The physical meaning of the unitary transformation (|38|) is that the number state of ¥(t)b(t) 
is a squeezed state of the corresponding one of a'(t)a(t). Similarly the number state of 

aJ(t)a(t) is a squeezed state of the corresponding one of a (t)a(t). 

Without losing generality, we shall use the quantum action operator 

lit) = ha\t)a(t) = hh(t), (44) 

which corresponds to the action variable (|9]). By writing the complex solution to Eq. (|5J) as 



(45) 



Eq. (H) becomes a nonlinear auxiliary equation 



c + 



XZ-Y 2 + 



XY -XY 



X 



+ 



X 3X 5 



2X AX 2 



c = l 



(46) 



where we used 0\ = l/( 2 from Eq. (|i~7|) . Then the invariant operator (p) takes the form 

m 
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Another form of the complex solution, 



u{t) = -j=at)e- ie ^\ 



(48) 



with 0£ = X/£ , leads to the form JTT 



lit) 



(49) 



where £ satisfies the auxiliary equation 



X 



L-^+[XZ-Y 2 + 



XY-XY\^ X 2 



X 



£ 3 



(50) 



The invariant operator (44]) is nothing but a number operator. Each eigenstate of Af is a 
number state 



n(t)\n, t) = n\n, t). 

The wave function for the number state defined by Eq. 
phase factor, is (see Appendix B) 



(51) 

lj), up to an arbitrary constant 



1 



1/4 



u 



(2n+l)/2 



(2h) n n\ 



2nhu*u 



\Ju*u. 



n \V2h 



u u , 



exp 



u 



2hX \ u* 



-yW 



(52) 



where H n is the Hermite polynomial. It satisfies the time-dependent Schrodinger equation. 
In the sense that one finds the Fock space of exact quantum states of the Schrodinger 
equation, even the time-dependent quantum oscillator (|l]) is exactly solvable in terms of the 
complex classical solution u(t). 

Finally we turn to the phase operator conjugate to the action operator (|44]). A phase 
operator should satisfy the commutation relation 

[i(t),B{t)]=ih, (53) 

where 9 is restricted to (0, 2ir). The commutation relation leads to the Lerner criterion 

[e i§{t \l(t)]=he i§(t) . (54) 

There are many schemes to define the phase operator, and each method has its own advan- 
tage and disadvantage 0. As there is a formal similarity between time-independent and 
time-dependent quantum oscillators, we follow straightforwardly the various well-known def- 
initions of phase operator. First, a naive definition from ( pT|) according to the correspondence 
principle would be 
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a(t) = e*°®nV 2 (t), tf(t) = n 1 / 2 (t)e^ D( * ) , (55) 



which is the procedure taken by Dirac ||. The #d formally satisfies the Lerner criterion (|5~T|). 
but only if one ignores the fact that to get Eq. ([54]) formally, one must multiply Eq. fl35|) 
by the singular operator n _1//2 . Furthermore, #d(£) is not a Hermitian operator, so e l0D ^ is 
not unitary 0. Second, the phase operator by Susskind and Glogower [E| is defined as 

E(t) = [n(t) + ir 1/2 a(t), E\t) = a\t)[h(t) + 1]~ 1/2 . (56) 

E and E^ are the analogs of e 10 and e~ 10 , which satisfy the Lerner criterion The 

analogs of cos# and sin# can be defined similarly. E and E^ have the equivalent number 
state representation 

oo oo 

E(t) = Y l \n,t)(n + l,t\, E\t) = £ \n+ l,t)(n,t\. (57) 

n=0 n=0 

However, E is a one-sided unitary operator, 

E(t)E\t) = l, E\t)E{t) = 1 - |0,t)(0,t|. (58) 

This is somewhat related with the ambiguity in defining a phase for a(t) = in classical 
theory or for the vacuum state with a(t)\0,t) = in quantum theory. Third, the phase 
operator defined by Pegg and Barnett in a finite dimensional subspace of Fock space |13| , is 
now given by 

s 

<W*)= J2 e m\0m,t)(e m ,t\, (59) 
m=0 



where 



\e m ,t) = ^L=ire m0 -\n,t), 

v s + 1 n=0 



2/TZ 7T 

9m = 9q ~\ r, (^o = constant). (60) 

s ~\~ X 



Fourth, the phase operator corresponding to the classical phase variable ( f22|) is naively given 
by 

^(t) = ^[lnat(t)-lna(t)]. (61) 

The phase operator ([61] ) formally satisfies the commutation relation (|53D and is a general- 
ization of the time-independent one in Ref. [0 to the time-dependent one. However, it is 
not really well-defined, since a has a zero eigenstate, so that its logarithm is a divergent 
operator. 

A few comments are in order. At present there is no consistent definition of a phase 
operator for a quantum oscillator free from all conceptual and technical problems 0. This 
dilemma is somewhat rooted in the quantization of phase variable according to Eq. (|53|) : 
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J = z^. (62) 

The eigenfunctions of Eq. (|62"D, consistent with the periodicity of 6, are given by 

V n (0) = -=e- in6 , (n = all integers). (63) 

The wave functions (|63"D form a basis of the Hilbert space on a circle and make the closure 
complete when n runs for all integers. Denoting the wave functions (|5"3"D by \n, 6), the 



following phase operator [14 



e i§ = ]T \n,9)(n+l,9\, e~ i§ = >T |n + 1, 0)(ra, 0| (64) 

n=— oo n=— oo 

becomes unitary 

e i§ e~ i§ = e~ i§ e i§ = 1. (65) 



However, compatibility with Eqs. fl5"T|) restricts n to all non-negative integers. As Eq. ( |52| ) 
and ( |63"D are the coordinate and phase- representations of the action eigenstate fl5TP, the 
phase operator (|>|) is the same obtained by extending Eq. ([57]) over all integers including 
negative ones. But there is no physical meaning to the wave functions of oscillator with 
negative eigenvalues. 



IV. CONCLUSION 

In summary, we found pairs of real and complex invariants for a time-dependent classical 
oscillator, which are linear in position and momentum and transform canonically the phase 
space to invariant spaces. An action was defined that is quadratic in each invariant space, and 
its conjugate phase (angle) variable was found. The relation among different phase variables 
was exploited. Analogously, pairs of linear complex invariant operators were found for the 
time-dependent quantum oscillator. The action operator, which is the quantum analog 
of the classical action, may be used to construct a Fock space of exact quantum number 
eigenstates. As the invariant operators act as the annihilation and creation operators, we 
followed the procedure for time-independent oscillator to define the phase operator conjugate 
to the action operator for the time-dependent oscillator. 

The action-phase variables were defined for the time-dependent classical oscillator exactly 
as for the time-independent oscillator. However any known scheme to define a consistent 
phase operator, even for the time-independent oscillator, confronts some technical difficulty 
@. The same degree of difficulty remains in defining the phase operator for the time- 
dependent oscillator. In this paper we did not intend to resolve all the puzzling issues of 
phase operator but to extend the current methods to the time-dependent oscillator. One 
interesting possibility will be to make use of the Wigner function of the exact quantum state. 
The Wigner function can be expressed in the invariant space, whose integral along each ray 
of constant phase may lead to a phase pseudo-probability distribution [I5]|. Then the phase 



probability distribution can be used to define an ensemble average of any physical function 
of the phase variable in quantum theory. 
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APPENDIX A: TWO LINEAR INVARIANTS 



In this appendix we derive two linear invariants given by Eq. or (|T6|). The Hamilto- 
nian (|I|) can be written as 

H(t) = V + (Al) 
where P and Q are new canonical variables 

« = 7Jft* P ^i" + W) q \ (A2) 



and 

u\t) =X(t)Z(t) -Y 2 (t). (A3) 
The Q and P are canonical conjugates of each other, satisfying 

{Q,P} PB = 1. (A4) 
By introducing another pair of new canonical variables a\(t) and a^it), satisfying 

{a 1 (t),a 2 (t)}p B = h (A5) 
we can express Q and P in terms of a\(t) and a 2 (t) 

Q = — [u 2 (t)ax(t) - u 1 (t)a 2 it)], 

P = -^=[u 2 (t)ai(t) - ni(t)a 2 (t)], (A6) 
X(t) 



and vice versa 



fx(tj JW) 



a 2 (t) = ^Lp-^L Q . (A7) 



Then the Hamiltonian has the form 
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H(t) = -2H Q (t)a x {t)a 2 {t) + #i(t)a?(t) + H 2 (t)a 2 2 (t), 



(A8) 



where 



^o(t) 
Hi(t) 



2X(t) 



\ui(t)ii 2 (t) + u 2 (t) Ul (t)u 2 (t)l 
[ul(t)+uj 2 (t)ul(t)}, 



2X(t) 1 

Finally we require ai(t) and a 2 (t) each to satisfy the invariant equation (@): 
= JUi(*) + {ai(f),#(*)}p B 



-{s(l) + r- y " + 

2 , XY-XY\(u x 



+ 



X 



X 



dl(t) 

a 2 (t), 



= T-a 2 (t) + {a 2 (t) } H(t)} PB 
at 



Uy 



d (ii 2 
Jt\X / 



XZ -Y z + 



, , XY -XY\ (u 2 



X 



X 



+ 



ay[t) 
a 2 (t). 



(A9) 



(A10) 



Therefore, a x (t) and a 2 {t) satisfy Eq. (|J) provided that and u 2 {t) satisfy the classical 
equation of motion (pj). 



APPENDIX B: WAVE FUNCTIONS 



In this appendix we derive the exact wave functions (|5^) for the time-dependent oscillator 
([[]) that satisfy the time-dependent Schrodinger equation 



ih^(q,t)=H(t)H>(q,t). 



(Bl) 



First we note that the operators a(t) and a^(t) in Eq. (|34l) are time-dependent analogs of 
time-independent annihilation and creation operators. In fact, in the time- independent case 
with constant X Q , Y , and Z , by choosing the complex solution to Eq. (|5|) given by 



u{t) 



-iijQt 



w o — X Z — Yq, 



(B2) 



it can be shown that d(t) = e lw °'do and a^(t) = e lulot aQ, where do and are time- 
independent annihilation and creation operators. The time-dependent phase factors are 
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necessary to satisfy Eq. ( pi|) . Hence n(t) = no and number states are the same in both 
representations. We follow exactly the same procedure of the time-indepenent oscillator to 
obtain the wave functions for the time-dependent oscillator. The ground state wave function 
is annihilated by a(t) 



(B3) 



so we obtain 



i 



1/4 



U 



1/2 



2irhu*u 



\Ju*u 



cxp 



u 



2hX V u* 



Y W 



(B4) 



Here we chose the time- dependent phase factor (u/y/vFu) 1 ^ 2 so that Eq. ( p4|) satisfies Eq. 
(|BT|). The nth wave function for the nth excited state is obtained by applying a*(t) n times: 



* n ( q ,t) = ^(tf(t)\ Mq,t) 



i 



1/4 



U 



(2n+l)/2 



(2h) n n\ 



2nTiu*u 



y/vF 



H, 



u . 



n \V2h 



exp 



u u . 



2hX\u 



h-rW 



■ (B5) 



It should be remarked that the wave function (|B5| ) satisfies not only Eq. (|B"l| ) but also the 
relation 



dqK(q,t)U^)y n (q,t) = I dq** n (q,t)H(t)* n (q,t). 



(B6) 



Hence Eq. ( B6|) implies that the geometric phase factor equals to the dynamical phase 
factor. In the special case of Y — 0, the wave function ( [Bo]) reduces, up to the phase factor 
(u/Vu*u) 1/2 , to Eq. (3.9) of Ref. §. 
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